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Q\ • We consider the flavour non-singlet Reggeon within the context of perturbative QCD. This 
consists of ladders built out of "reggeized" quarks. We propose a method for the numerical 

i-^ I solution of the integro-differential equation for the amplitude describing the exchange of 

such a Reggeon. The solution is known to have a sharp rise at low values of Bjorken-x when 
applied to non-singlet quantities in deep-inelastic scattering. We show that when the running 

^ ■ of the coupling is taken into account this sharp rise is further enhanced, although the Q^ 

dependence is suppressed by the introduction of the running coupling. We also investigate 

k> i the effects of simulating non-perturbative physics by introducing a constituent mass for the 

^ ' soft quarks and an effective mass for the soft gluons exchanged in the t-channel. 
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1 Introduction 

In recent years, much attention has been given to the perturbative QCD simulation of the 
Pomeron. The reason for this is that the reach of HERA is such that one now has data on 
structure functions and differential cross-sections for other inclusive processes which are well 
into the diffractive region (low-x) whilst at the same time maintaining momentum scales for 
all the kinematic variables. These variables are large enough that a renormalization-group 
improved perturbative expansion summed to all orders in leading Inx is expected to be 
valid. This is the region in which one expects to be able to test the BFKL Pomeron [|l|. 
Notwithstanding this, it should be recalled that the original motivation for the study of the 
Pomeron in QCD was an attempt to explain how the successes of Regge theory could be 
under-written by a quantum field theory. 

Somewhat less emphasis has been placed on the Regge trajectory below the Pomeron (the 
Reggeon) for which phenomenological fits [0 have shown have an intercept ~ ^. In the 
same way that the QCD Pomeron is constructed from ladders of reggeized gluons with a 
colour singlet projection, so the QCD Reggeon is constructed from ladders of a reggeized 
quark-antiquark pair, again with colour singlet projection. 

The cleanest way to distinguish between Pomeron dominated processes and Reggeon dom- 
inated processes is to note that a quark-antiquark pair can be in a flavour non-singlet 
state. We therefore consider quantities which are controlled by flavour non-singlet oper- 
ators. Such quantities include the structure function F3 in deep-inelastic scattering, or the 
spin-dependent structure functions. The Reggeon is only expected to dominate such quanti- 
ties at sufficiently low-x, and the extraction of low-x data for such quantities is very difficult. 
The spin-dependent structure functions at low-x have recently been considered in ref . [Q . At 
the moment the lowest values of x for which we have data on F3 are not yet in the asymptotic 
region where we would expect the Reggeon to dominate, but on the other hand they are not 
very far away and in the same way that HERA is currently used to test the QCD Pomeron, 
the QCD Reggeon could reasonably be expected to be probed in the not-too-distant future. 

In order to construct the Reggeon within the context of perturbative QCD, one needs first to 
establish that the quark reggeizes, in the same way that the gluon reggeizes. In other words, 
it is necessary to show that to leading order in Ins, the amplitude in which the quantum 
numbers of a quark is exchanged in the t-channel, has an s dependence given by 

From standard partial wave analysis we expect aQ{t) = | in leading order and for the t 
dependence to arise from higher order corrections. This was first established some time ago 
0. However, we feel that several important features of this derivation were not sufficiently 
clearly explained in ref. ||^ . We have repeated the relevant calculations and we briefly review 
our method for deriving the reggeization of the quark in section 2. Once quark reggeization 
has been established, an integro-differential equation can be established for the sum of the 
leading In s parts of the amplitude for the exchange of the quantum numbers of the Reggeon. 
This has been derived in ref. and we quote the result in section 3. This equation can, in 
principle, be solved analytically since it involves a kernel which is conformally invariant (in 
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Figure 2.1: The tree- level exchange of a soft quark 

two dimensions) such that its eigenfunctions are the representations of the conformal group. 
Unlike the case of the Pomeron the leading eigenvalue is not analytic in the coupling. This 
causes severe difficulties when one tries to reproduce the analytic results by numerical means. 
We propose a programme for numerical solution of the integro-differential equation which is 
consistent with the known asymptotic behaviour. This allows one to introduce modifications 
to the kernel, such as the running of the coupling, which breaks conformal invariance, so 
that the analytic approach is no longer viable, but the numerical approach remains reliable. 

In section 4 we apply this to a model of fiavour non-singlet deep-inelastic scattering. We 
call this a model since an arbitrary function has to be taken for the impact factor describing 
the emission of a Reggeon from the target hadron. However, we do not expect either the x 
dependence at sufficiently low x, or the Q^ at sufficiently high Q^ to be particularly sensitive 
to the exact nature of this impact factor. We look at the effect of running the coupling and 
observe that whereas this further enhances the predicted rise in the non-singlet structure 
functions as x — ;> 0, it suppresses the Q^ dependence of such structure functions. Finally, 
we postulate that the most significant feature of non-perturbative effects in the case of 
the Reggeon is to provide the soft quark and antiquark exchanged in the t-channel with a 
constituent mass, and the soft gluons an effective mass. We therefore solve the modified 
equation for the Reggeon in which a mass is assigned to the soft particles. We find that 
this significantly suppresses the low- a; rise and that this suppression extends to values of 
Q^ way beyond the assigned values of the constituent mass. In section 5 we summarize our 
conclusions. 



2 The Reggeized Quark 

In this section we review the derivation of the "reggeized quark" . The result obtained agrees 
with that previously obtained [^], but in our derivation some of the subtleties leading to the 
result are discussed in more detail. 

We begin by considering the process 

q + q -^ g + g, (2.1) 

in the Regge regime, s ^ |t|, where a soft quark is exchanged in the t-channel as shown in 
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Figure 2.2: One-loop leading Ins corrections toq+q—>-g + g. The 2-body cuts in the 
s-channel are shown, however (iii) and (iv) have similar cuts in the w-channel 



Fig. prri. The amplitude is given by 



A(g) 
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(2.2) 



where we have introduced holomorphic coordinates in the plane transverse to the incident 
momenta pi and p2 such that 7 (7*) are the 7-matrices in this plane, q (q*) are the components 
of the transverse momentum transferred in the t-channel, and ei (e*) and 62 (£2) ^ire the 
polarisation vectors of the outgoing gluons. 

We now consider the one-loop corrections which are proportional to Ins, since it is only 
these contributions that are relevant for reggeization. The relevant diagrams are shown in 
Fig. p.2| , and each of these gives the same contribution (up to a colour factor), which may 
be written as Q 
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(2.3) 



The diagrams of Fig. p.2| have colour factors 



^ The quantity k^ , used to scale s inside the logarithms is understood to be a square momentum which is 
of the order of the square momentum transfer |i|, or some other momentum (^ s) involved in the coupling of 
the Reggeon at the top or bottom of the ladder. Since we are confining our discussion to leading logarithm, 
its exact value is unimportant. 
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Figure 2.3: The effective quarlc-gluon vertex 
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where A^ (= 3) is the number of colours. 

Adding these one-loop contributions we obtain the result for the one-loop correction to be 

A(g)=A(g)ln(^)eQ(g) (2.4) 

with 



eQ(g) 



-C, 
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(2.5) 



and Cf = {N^ - l)/2N. 



At the next order a number of complications occur, some of which are common to the problem 
of gluon reggeization ||T| . The first of these is that the diagrams that contribute in leading 
In s are not all of the ladder type. However, all diagrams whose s-channel cuts contain three 
intermediate particles can be reduced to a ladder-type diagram using an effective vertex. 



This is demonstrated in Fig. p.3| . The Feynman rule for this effective vertex involves the 
substitution of 7°^ in the standard rule with Fg, where 
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from which it follows that 
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+ h.c. 



(2.6) 



(2.7) 



rr(fc., h-i) ■ Q = -^ ( ' 7 + h.c. ) . {21 



/Cj_i, ki are the transverse momenta of the incoming and outgoing quarks, defined in terms 
of Sudakov variables by 

with k\i = —\ki\^. 

We also need the equivalent effective vertex for the emission of a gluon with polarisation e^ 
from incoming and outgoing vertical gluon lines with momenta fcj-i, ki and Lorentz indices 
fj,, p. This is given by 

2p2Pi / ^iki-ik': 

In both of these effective vertices, we have assumed the multi-Regge kinematics, namely 

ai-i > ai 

m > iA-ii. 

The product of these two effective vertices, summed over the polarisations for the emitted 
gluon is given by 

V-/T. N/T. ^ f . k'*\q~k\'^ k*\q-k'\'^\ , , , 

E (rc^ ■ ^)(rQ ■^)=lU- i^'l^J - J_feJ J + h.c. (2.9) 

pol'ns ^ ^ 

A further complication arises from the need to consider diagrams in which the central cut line 
involves a quark or antiquark (as opposed to a gluon) and this is demonstrated pictorially 
in Fig. 



Piecing all this together with the appropriate colour factors we find that the contribution 
from all such diagrams may be written as 

Mq) = \{^Q{q)fMq) + Mq)' (2.10) 

where 



A2{q)' = -{c, + ^) C,{r\%,^^\n' {{r.)^{P2)lu{p,)elel 
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rf^kirf^ka— -^^ -—:t- p-p- + h.c. 

l^2rlQ' — ki\'^\ki — k2\ 



The first term suggests that the amplitude is indeed "reggeizing" , whereas the contribution 
^2(9)' cancels against the diagrams shown in Fig. p.5| , in which there are two particles in each 
(s-channel or w-channel) cut and where the loop on either side of this cut is interpreted as the 
reggeization of the vertical quark or gluon line. We have now developed a general mechanism 
for the derivation of the reggeized quark, using the traditional "bootstrap" approach. We 
make the ansatz that the quark reggeizes, with Regge trajectory given by 1/2 + eQ^q) as 
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Figure 2.4: Examples of 2-loop ladders. Diagrams with a quark or antiquark rung have to 
be considered in addition to the gluon rungs with effective vertices 



defined in ( |2.5| ) and show that this ansatz is self-consistent (the summand of 1/2 will be 
explained shortly). We also need to make use of the reggeization of the gluon which has a 
Regge trajectory of 1 + eg(g) with 



ee(g) 



2 27r2 



d^k 
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(2.11) 



In order to develop an expression for the amplitude for the process ( ^.11) , whose iterative 
solution gives a series 

Ao{q) + Ai{q) + A2{q) + ■ ■ ■ , 

we need to consider the addition of both gluon and quark or antiquark rungs. This is 
demonstrated in Fig. P]B. In general these ladder diagrams lead to an integro-differential 
equation of the form 



dF{s, k,q) 
dins 



f d^k']Cik,k\q)F{s,k',q), 



(2.12) 



in which the kernel, /C, may be divided into two parts. The effect of adding a rung gives a 
contribution of 
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Figure 2.5: The complete set of leading Ins two-loop, two particle cut diagrams 
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Figure 2.6: The addition of a rung involves the consideration of both a gluon rung with 
effective vertices and of a quark rung. Only the upper-half of the diagrams is shown, and 
note that the quark rung diagrams imply that we must have $1 = $2 



The first term is precisely what one would expect from the reggeization of the quark, whereas 
the second term cancels when we take into account the fact that the quark and gluon 
exchanged in the t-channel must themselves be reggeized and this reggeization gives rise to 
a correction 



j £\^lC^^^\k,k\q) = {eQ{k) + eg{q - k)) 



(2.14) 



When all is pieced together we find that the effects of adding a rung and of taking into 
account the reggeization of the quarks and gluons exchanged in the t-channel is such that 



d'k' (/C-'^s(fc, k', q) + IC^^^^^k, k', q)) = eQ{q), 



(2.15) 



which means that the leading In s part of the amplitude, summed to all orders in perturbation 
theory, gives rise to an amplitude of the form 



^e^'Uo 



~ s 



l/2+eQ(g) 



(2.16) 



where the "naive" s dependence s^/^ arises from the usual s dependence of an amplitude in 
which a spin-| particle is exchanged in the t-channel and is contained in the normalisation 



of the spinors in ( 2.2 



3 The Integral Equation for the Reggeon Amphtude 
and its Solution 

The integral equation for the perturbative amphtude for processes involving the exchange 
of a Reggeon are obtained by considering ladders in which the vertical lines are a reggeized 
quark-antiquark pair, and the horizontal rungs are gluons which couple to these reggeized 
quarks via the effective vertex given in (|2.7| ). The colour singlet is projected from the two 
fermions exchanged in the t-channel. This has been studied in [^. One obtains an integral 
equation for the Mellin transform of the amplitude, f{s, k, q) defined by 

/(w, k,q)= — \^—j f{s, k, g), 

where k is the transverse momentum of the incoming quark and q is the momentum transfer 
(which may be taken to be transverse). In terms of this Mellin transform one obtains 



"^ f ^2, / f k' ■ (k' - q) 



f(uj,k,q) = foiu;,k,q) + ^ f d'k' \ 



k'^(k'-q)2 



k k' 



k'k (k - q) ■ (k' - q) 1 ~, , , , 

'^ ^c^-k-|,,,^^+ i^-l-f-l K(^.M), (3.1) 



27rcjy ik'2(k-k')2 (k' - q)2(k - k')2 

where we have used the notation 

as = Cp — , 

IT 

and foiuj, k, q) is the Born term obtained from the leading order contribution consisting of 
the exchange of a quark-antiquark pair (with no gluons). 

There is an important difference between this equation and the equivalent equation [|l| for 
the Pomeron, which is built out of reggeized gluons in a similar fashion. In the case of the 
Pomeron a solution which is of the reggeized form (i.e. a simple pole in the Mellin transform 
at a value of u equal to the Regge trajectory) is prevented by the fact that the colour factor 
which arises in the projection of a colour singlet from the two gluons in the t-channel is 
twice that obtained from the projection of a colour octet (relevant for the integral equation 
describing the reggeized gluon itself). In the case of the Reggeon, a colour factor of Cp is 
obtained both in the case of the reggeized gluon, and in the projection of a colour singlet 
from a quark-antiquark pair. In the Reggeon case, simple reggeization is prevented by the 
uj dependent factor 

k k'^ 



which occurs inside the kernel on the RHS of ( pj.l| ). This term arises from the fact that when 
one exchanges fermions in the i-channel such that the momenta of these fermions appear 
in the numerators of the propagators, care must be taken in setting the kinematic limits 
of the phase-space integral over the longitudinal components of momentum p[. It is these 
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kinematic limits which lead to the uo dependent term in the kernel and, as we shall see later, 
this has a radical effect on the solution to this equation. 

Henceforth we shall be confining our discussion to the case of zero momentum transfer 
(g = 0), for which for which the second term on the RHS of ( |3.1|) simplifies to 



ICiuj, k, k') (g) f{u, k') 



a„ 



dk" 



mm 



k k' 

k^ + k'^ 



/>,fc') 
f{uj,k')-~f{uj,k)]+f{uj,k)\, (3.2) 



where we have integrated over the angular part of /c, assuming that the Reggeon is dominated 
by an amplitude which is azimuthally symmetric, i.e. that f{uj,k') depends only on the 
magnitude, \k'\. 

This can be solved exactly by exploiting the two-dimensional conformal invariance of the 
kernel to conclude that a function of the form 

f{u,k) = {\k\T 

is an eigenfunction of the kernel with eigenvalue ujIu) given by the solution to the equation 



a. 



ijj 



u 



- 2^(1 + iu) - 2^(1 - lu) - A-fi 



(3.3) 



_cuV4 + z/2 
The leading eigenvalue is given by 

ujQ = a/25^, 
which implies that the Reggeon has an asymptotic s behaviour of the form 

Another qualitative difference between the QCD Reggeon and the QCD Pomeron lies in 
the fact that such a leading s behaviour, which is not analytic in the coupling, cannot be 
reproduced in an ordinary perturbative expansion in a^, although the even orders in the 
expansion of this s dependence matches the double leading logarithms found in ordinary 
perturbation theory |^. 

Thus, in principle, we have an analytic method for solving the integral equation for the 
Mellin transform of the Reggeon exchange amplitude. However, in order to be able to 
introduce further refinements, such as the running of the coupling or some simulation of 
non-perturbative effects which necessarily break the conformal invariance, one needs to be 



able to carry out a programme of numerical solution of (|3TT 



^ We note here that — wq is also an eigenvalue of the kernel, and although it is sub-leading, it will turn 
out to be convenient to exploit this in our subsequent analysis. 
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This is most easily achieved by inverting the Melhn transform to obtain an integro-differential 
equation 

^^''^'^^ -)C{s,k,k')(g)f{s,k'), (3.4) 



91ns 



where 



"1 dz 1z^ 
lC{s, k, k') ® f{s, k') =&s{ j -Y—i (^(^' ^^) + ^(^' ^/^) - 2/(3, fc)) 

-{j{szM)^f{sz,klz))\ (3.5) 

l/s ^ J 

and we have introduced the variable z to represent min(fc/fc', k' /k). The uj dependent part 
of the kernel in Mellin transform space is now encoded by the scaling of s by a factor of z in 
the appropriate terms. 

We know that the asymptotic behaviour of the solution to this differential equation should be 
s'^°. The usual method of numerical solution of such differential equations is to start with an 
initial function which has no s dependence and to integrate in steps using the Runge-Kutta 
method. Unfortunately, since the leading s behaviour arises from the u dependent part of the 
kernel in Mellin transform space, which has now been transferred to the s dependent part of 
/(s, k) on the RHS of (|3.5| ), an initial function which does not contain any s dependence can 
never generate the expected asymptotic behaviour in any step-by-step numerical integration 
routine. 

This difficulty can be circumvented by extracting the leading s'^° behaviour. This is achieved 
by defining a function g{s,k) by 

f{s, k) = cosh (tuo ln(s)) g{s, k). (3.6) 

The asymptotic behaviour of /(s, k) is now displayed explicitly and we expect the function 
g{s, k) to have a far less dramatic s dependence. The integro-differential equation for g{s, k) 
is given by 

f/fl ( S n I 

— = sech{uQ\n s)IC{s, k, k') ® {cosh{uJo\n s)g{s,k')) — uJotai]h{LUQ\ns)g{s,k) (3.7) 
a Ins 

where the kernel )C{s,k,k') is defined in ( |3.5| ). We note here that g{s,k) = constant is an 
exact solution of (|3.7|) . We may therefore start with a function, go{k), that depends only on 
k but not on s (which one would normally obtain from an "impact factor" which encodes 
the coupling of the Reggeon to the particles between which the Reggeon is exchanged) and 
generate (after sufficient iterations) a solution for the Reggeon exchange amplitude with 
initial value go{k) taken as the value of the amplitude at low s. 

We have checked numerically that this method reproduces the known analytic results in the 
case of a trial initial function, which has a simple Fourier transform. We have taken 

^o(fc)= ,i;i 1^2 =- / duT{l+zu)T{l-zu){\k\Y. (3.8) 
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Figure 4.1: The impact factors $^ and $p encode how the virtual photon (or W-, Z-boson) 
and the proton respectively couple to the Reggeized quark ladder. The central circle sym- 
bolises the kernel's role in adding rungs and Reggeizing the uprights 



For large s the amplitude is then given by 



1 



f{s,k) = - / duT{l+iu)T{l-iu)(\k\^y''s'^^'\ 



(3.9) 



where (^{u) is given by ( |3.3| ). This can be compared with the numerical method described 
above. Having established the validity of this numerical approach, one can introduce re- 
finements to the kernel, which render it intractable to analytical methods, but still open to 
numerical analysis. In the next section we discuss the results of such an approach. 



4 Numerical Results 

In this section we apply the technique discussed in the last section to a model of deep- 
inelastic scattering processes at low-x, which would be dominated by Reggeon rather than 
Pomeron exchange (such as the structure function F^ for virtual W- or Z-boson scattering, 
or the spin-dependent structure functions which are intrinsically flavour non-singlet). For 
the rest of this section we shall use the term "structure function" to mean the flavour non- 
singlet part of a structure function, which we expect to be dominated by Reggeon exchange 
at low-x. In order to do this we need some model for the "impact factors" which account for 



the coupling of the Reggeon at the top and bottom of the ladder, as indicated in Fig. 4.1 



Unlike the case of the Pomeron, there is a direct coupling of the quarks in the Reggeon to 
the virtual particle involved in the scattering (photon, PF-boson, or Z-boson). This means 
that the upper impact factor is proportional to 
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Figure 4.2: The (lack of) dependence of results on the parameter rj in ([4.1|) is demonstrated. 
The structure function at fixed Q^ is given for six values of 7] from 0.6 to 2.0 



where —Q^ is the square momentum of the virtual boson. At the other end of the ladder 
we need the impact factor for the coupling of the Reggeon to the target proton. This is 
non-perturbative in nature and so we model this with the simple form 



$p oc e 



-k^/rjA^ 



QCD 



(4.1) 



where we expect r] to be of order unity, and we have not fixed the constant of proportionality 
which means that our scale remains arbitrary. This is clearly a crude estimate of the lower 
impact factor, but we do not expect our results for the low-x, large Q^ regime to depend 
critically on the exact form of this impact factor. This is because after a sufficient number 
of iterations of the kernel one expects the amplitude to lose all "memory" of the initial 
condition. As a demonstration of this, we show in Fig. [4.2| the variation of the structure 
function plotted against x for Q^ ~ 100 GeV^ for a range of the parameter t] between 0.6 
and 2.0. We see that the variation of the shape of the structure is indeed small and the 
differences can be substantially absorbed into an overall normalization. 

We start by assuming fixed coupling. In this case the behaviour of the structure function at 
low X for each value of Q^ was obtained by starting with an initial function given by ( [4.1|) 
and numerically integrating ( ^.7] ) using the Runge-Kutta method with the coupling a^ taken 
to be as(Q^). From the resulting function /(x, k) (s is replaced by l/x), the point k'^ = Q^ 
was selected, thereby imposing the delta-function which represents the impact factor at the 
top of the ladder. This process was then repeated for different values of Q^. 

Next we wish to include the effect of the running of the coupling. This is achieved by 
promoting the coupling a^ in the expression for the kernel to a running coupling. We take 
the larger of the momenta k and k' as the argument of the running of the coupling, which 
means that the coupling must now be taken inside the integral in (|3.5|) . One immediate 
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effect of this is that the quantity Uq, which was the leading eigenvalue of the kernel in the 
case of fixed coupling, now becomes a function of k given by 



/2 



1 f dk'' 
^o{k) = - — 2-as(max(A;, /c')) 



k k' 



(4.2) 



and it is this value of uj{k) that is used to factor off the leading behaviour. Thus we define 
g{x, k) as 

/(x, k) = cosh {uj{k) ln(s)) g{x, k), 

and the x dependence of g{x, k) is given by ( p.7|) with cxg replaced by a;s(max(A;, k')). Note 
that in this case, unlike the case of fixed coupling, we only need to solve the corresponding 
differential equation once, since the argument of the coupling for each iteration does not 
now depend on Q^. A single pass through the range of x returns a function /(x, k) and the 
results for the structure function at any value of Q^ can be read off by setting k"^ = Q"^. 

As is always the case when we introduce a running coupling in such integro-differential 
equations, we need to impose an infra-red cutoff below which a^ takes the value Omax and 
ceases to run. This is because the integral in (|3.5|) samples all possible momenta so that 
without the imposition of such an infra-red cutoff one would run into the Landau pole. We 
have chosen am»^ to take the value 1. 



In Figs. |1]^ and ^^ , we show the resulting structure function /(x, Q ) plotted against x for 
different values of Q^. Fig. |4.3| refers to fixed coupling and Fig. ^^ to running coupling. 
We have also displayed a curve indicating the x~^/^ dependence that one would expect to 
obtain from the naive counting of amplitudes in which spin-| particles are exchanged in the 
t-channel. The enhancement at low-x is clearly seen and at sufficiently low-x the structure 
function behaves like 

for the case of fixed coupling, and for running coupling we have an even further enhancement 
of the low-x rise which asymtotically behaves like 



On the other hand, we see from Figs. ^^ and ^^ that there is much less Q dependence in 



the rise at low- a; when the running of the coupling is taken into account. This is consistent 
with recent results obtained in the case of the Pomeron with running coupling. This can 
be attributed to the fact that for the Q^ dependence of the structure functions at low-x, the 
effect of running the coupling is to sample values of k'^ which are larger than Q^ for which the 
running coupling is smaller than the fixed coupling value, whereas for the Q^ independent 
part one is also probing the infrared region where the running coupling is enhanced. 

Finally, we consider possible non-perturbative effects. One important consequence of non- 
perturbative QCD is that quarks acquire a constituent mass. As soft quarks are exchanged 
in the t-channel in the Reggeon, it is quite likely that this will be the most significant 
contribution of non-perturbative QCD to the behaviour of the Reggeon. We also choose 
to give the soft gluons an effective mass, and so to simulate the effects of non-perturbative 
QCD we assigning a mass, m, to the soft quarks and gluons. 
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Figure 4.3: The rise in the structure function at low-x is shown here for four different values 
of Q (in [GeV]) for static as- The naive dependence of x"!/^ is also shown for comparison 
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Figure 4.4: The same as figure ^^ but with a running coupling. Note that the low- a; rise is 
now even more pronounced 
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Figure 4.5: The structure function / shown against Q for three values of x. The couphng 
has been fixed for this calculation 







40 50 

Q [GeV] 



80 90 



Figure 4.6: Similar to figure [4.5| , but with the coupling now running. This acts to reduce 
the rise at small Q 
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The kernel fB.5|) then becomes 



)CUs,k,k')(^fis,k') 



^ dz 



as{k'^ 



1+z'^ + M^ 



z^ + M^ y ^((1 + zy + M2)((i - zy + m^) 

Xjz" 



X 



1 + 1/^2 ^ ^2 



+ 



v/((l + \lzf + M2)((l - Ijzf + M2) 



1/^2 + ]Vf2 

l|(/(.,A;/^)-/(s,fc)) 



l/s 



as(/c 



Z2 + M2 



/(sz,A;z) + a,(A;7z^) 



1/^2 



1/^2 + ]Vf2 



f{sz,k/z) 



(4.3) 



in which we have defined M2 = rn^ /k"^. 
The expression for ct;o(A;) then becomes 



a;o(fc^ 



(i/;; 



/2 



2 ;./2^ 



2i A;'2 + m2 



Q;s(max(/c , /c' )) 



mm 



A; A;' 



i^oik^) 



(4.4) 



and it is this expression for ujQ{k'^) which is used to define g{s,k) in ( |3.6| ). 



The introduction of the mass term acts to reduce the value of ujq, and a plot of this effect 
is given in Fig. |4.7|. Note that the smaller Q^ values receive a more marked decrease in the 



leading eigenvalue, which further acts to limit growth in the infra-red region. 
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Figure 4.7: The variation of the leading eigenvalue with propagator mass, shown for four 
different values of Q^ (all given in GeV^) 
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Figure 4.8: The effect of non-zero propagator masses is shown at z = 10~^. Increasing the 
mass (all in GeV) leads to a reduction at large Q due to the decrease in the leading eigenvalue 

UJQ 
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Figure 4.9: The same data as in figure [4.8| is shown here with a logarithmic x-axis to 
demonstrate the infra-red effects of the propagator mass 
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We show the effect of the massive propagators in Fig. [4.8| . Here we have plotted against Q for 
five values of m, including ?7i = 0. The effect is quite strong, even at larger values of Q, with 
the reduction in / being more pronounced than might be expected of a measure originally 
designed to affect just the infra-red region. The reason behind this has of course just been 
detailed: increasing m reduces ujq so that by the time we have evolved the structure function 
to low-x (the graph is at a; = 10"^) it is substantially reduced across its whole domain due 
to the weaker leading behaviour. If we were to consider the same function at larger x the 
differences between the graphs for differing m would not be so pronounced for large Q^. 



In Fig. O we have shown exactly the same data as in the previous graph, but with the Q 
scale now logarithmic in order to illustrate the small Q region more clearly. As expected, 
the introduction of the mass term results in a reduction at low Q, but whereas the larger Q 
reduction only becomes apparent as we evolve to low-x, the growth in the Q ^ m region is 
immediately regulated. 







Figure 4.10: The upper set of four lines show zero mass results to provide comparison with 
the lower set which have m = 0.3 GeV. As expected the rise at low x has been slowed by 
the non-zero mass 



In Fig. |4.10| we have again plotted the zero mass/running coupling graphs for four Q values, 
but we have also included on the same plot the case where we take m = 0.3 GeV. Thus the 
upper set of lines are the same as in figure [4.4| , and the lower set show the relative effect of 
the propagator's mass. As would be expected from our previous arguments the low-x rise in 
/ is reduced. 
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5 Summary 

We have reviewed the derivation of the reggeization of the quark. This reggeized quark is 
used to construct the integro-differential equation for the amphtude of a process in which 
the quantum numbers of the Reggeon are exchanged. We have proposed a method for the 
numerical solution of this equation which is consistent with the leading large s behaviour. 
This enables one to make modifications to the kernel in order to account for a running 
coupling and for simulations of non-perturbative effects. 

The experimental quantities which can be used for probing the QCD simulated Reggeon are 
flavour non-singlet quantities at low-x. In all probability the reach of current experiments 
in deep-inelastic scattering are not yet sufficient provide us with such a probe, but it may 
well become possible in the near future. 

We find from our numerical studies that the introduction of the running of the coupling 
enhances the already sharp rise expected in the structure functions as x — >■ 0. On the other 
hand the Q^ dependence is considerably moderated when a running coupling is introduced. 
The effects of QCD beyond perturbation theory have been estimated by assigning a mass to 
the soft particles exchanged. We find that this has a dramatic effect in reducing the rise at 
low-x, even at values of Q^ which are considerably larger than the values of the constituent 
mass inserted. All of this suggests that it will be even harder to isolate the kinematic region 
where we might expect the results of perturbative QCD to dominate, and although a rise at 
low-x is still expected, it may well not be as dramatic as initially expected. 
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